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1. General framework 

Let t and x two real variables. We denote by IR^ the set of line vectors : 
(t, x) e IR^ and IR^'* = IR x* IR the set of column vectors that are the 
transposed from line vectors : 

(1.1) (*) e IR^.'. 

For e IR, we denote by the special Galileo transform defined by the 
relations 



(1-2) 



X I \ X — vt I \ —V 1 I \ X 



and we introduce also the space symmetry q defined by the conditions 

« (l) = i-x) = (0 -1) (!) ■ 



(1.3) 



Proposition 1. Relation between the generators. 

We have the following relations between the operators and q : 

(1.4) \fv,weM, Vv^Vw = Vy + w 

(1.5) q'^ = Id 

(1.6) G IR, yy •q^yy = q 

where Id denotes the identity matrix. 

Definition 1. Galileo group. 

We denote by GL^ (IR) the group of two by two invertible matrices with real 
coefficients. By definition, the Galileo group G is the subgroup of GL^{1R) 
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generated by the matrices yy for v G IR, the matrix q and the relations (1.4), 
(1.5) and (1.6). 

Definition 2. Space of states. 

Let m be a positive integer. The space of states Q, is an open convex cone 
included in IR"^' * : 

(1.7) f^ciR^'* and ^Wen, VA>0, xwen. 

Hypothesis 1. 

The Galileo group operates linearly on the space of states. 

For g E G and W E the action g^W of g on the state W is well defined 
as a linear operation of the group G on the set O : 

(1.8) ygeG, \fWen, g.Wefi 

(1.9) yg^g'eG, WW en, {g g')^W = g.{g' .W) . 

From an algebraic point of view, there exists an m by m set of matrices Y{v) 
and a matrix R with m lines and m columns such that 

(1.10) Vi;eiR, ywen, Y{v).wen 

(1.11) ywen, R.wen. 

Then the relations (1.4), (1.5) and (1.6) between the generators and the hypothe- 
ses (1.8) and (1.9) can be written in the vocabulary of the m by m matrices 

Y{v) and R : 

(1.12) Vt;, w;elR, Y{v) .Yiw) = Y{v + w) 

(1.13) = Id 

(1.14) V-uGlR, Y{v).R.Y{v) = R. 



2. Conservation laws 

We introduce a regular flux function Q 3 W i — y fiW) G IR"^'* and we 
consider the associated system of conservation laws in one space dimension : 

(2.1) g-W{t,x) + ^f{W{t,x)) = 

where the unknown function [0, +cx)[xlR 3 {t,x) i — y W{t,x) G Q takes 
its values inside the convex open cone Q. We suppose that the system (2.1) of 
conservation laws admits a mathematical entropy (Friedrichs and Lax [FL71]) 
r] : Q 3 W I — y ri{W) G H which is here supposed to be a regular (of class), 
strictly convex function admitting an associated entropy flux ^(•). Recall that 
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a pair of (mathematical entropy , entropy flux) allows any regular solution of 
system (2.1) to satisfy an additional conservation law (Godunov [G06I]) : 

(2.2) ^^rj{W{t,x)) + ^ [mit,^))] = 

and that a weak entropy solution of system (2.1) is by definition (see e.g. 
Godlewski-Raviart [GR96]) a weak solution that satisfies also the inequality 

(2.3) §-^vmt.,x)) + §^[mit,x))] < 

in the sense of distributions. 



• We introduce the Frechet derivative dri{W) of the entropy and its associated 
partial derivatives Q, 3 W 1 — y (p{W) G $ that takes its values in some set 
$ C IR"^ and defines the so-called entropy variables : 

(2.4) \fWen, Vr€lR"''\ dr]{W).r = (fi{W).r. 

We consider also the dual function of the entropy : ^ 3 (p 1 — y 77*((^) G IR in 
the sense proposed by Moreau [M066], that is 

(2.5) V*{^) = sup (ip.W - 7]{W)) . 

wen 

This dual function satisfies 

(2.6) dr]*{(p) = d(fW{(p) , where (f = dr]{W{(p)) . 
Hypothesis 2. Entropy flux and velocity. 

There exists a regular (of class) function u : Q 3 W 1 — > u{W) G IR which 
allows to write the entropy fiux ^(•) associated to the strictly convex entropy 
ri{») under the form 

(2.7) yWen, ^{W) = r]{W)u{W). 

• For the study of gas dynamics in Lagrangian coordinates, Despres [De98] 
has developed the general case where the function Q, 3 W 1 — > u{W) G IR 
is identically null. Hypothesis 2 is very little restrictive because we can define 
u{W) by the relation u{W) = as soon as the state W is such that 

the entropy r]{W) is not null. We observe also that field u{») is homogeneous 
to the ratio of space variable x over time variable t. For this reason, we have 
introduced in [Du99] the following 

Deflnition 3. Velocity of the state W. 

When hypothesis 2 is satisfied, we call velocity of the state W the scalar u{W) 
defined according to the relation (2.7). 
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Definition 4. Thermodynamic flux. 

When the hypothesis 2 is satisfied, the thermodynamic flux j : Q, 3 W i — > 
j{W) e IR"^'* is defined according to the relation 

(2.8) j{W) = f{W) - u{W)W. 

The hypothesis 2 is satisfied if and only if we have the following relation between 
the thermodynamic fiux and the velocity u{») : 

(2.9) VVFgO, <f{W).dj{W) + r]*{(f{W)) du{W) = 0. 
Proof of Proposition 2. 

It is sufficient to express the classical compatibility relation (see e.g. [GR96]) 
between the entropy fiux, the derivative of entropy and the derivative of the 
flux function, i.e. d(77u) = (p»df; this relation is what is usefull to deduce 
the relation (2.2) from the original system (2.1) of conservation laws. We derive 
the relation (2.8). Then we get df{W) = {du{W)) W + u{W) dW + dj{W) 
and 

(fdf- diiju) = 

= {du{W))ip.W + u{W)if.dW + ^•djiW) - {ipu{W) + r]du{W)) 
= du{W) {ip.W - ri{W)) + ifdj{W) because dW = Id(IR"^'t) 

= Ti*{ip)du{W) + (p»dj{W) 
and the desired result is established. □ 



Example. Gas dynamics (i). 

For the Euler equations of gas dynamics, we have m = 3, 

(2.10) O = |(p,g,e)^€lR^'*, p>0, e>|^}- 

Classically, the velocity u and the internal energy e are introduced thanks to 
the relations 

(2.11) q = pu , e = pe -\- - pu 

and because the pair (mathematical entropy , entropy flux) is of the form 

(2.12) (77,0 = i-ps, -psu) 

where the entropy s is the speciflc thermostatic entropy (see e.g. [Du90]), the 
relation (2.7) is satisfled and the deflnition 3 is perfectly compatible with this 
classical physical model for gas dynamics. Moreover, the classical relation be- 
tween the massic internal energy e , the massic entropy s and the speciflc volume 
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(2.13) de = Tds - pdr 

allows us to define the temperature T and the thermodynamic static pressure p 
(see e.g. Callen [Ca85]). Then the Gibbs-Duhem relation associated to extensive 
fields shows that the specific chemical potential fi can be defined as a function 
of temperature and of static pressure according to the relation 

(2.14) fi = ii{T, p) = e - T s + pT . 

After some lines of easy algebra, it is possible to express the entropy variables 
for the system of Euler equations for gas dynamics : 

(2.15) = i(^ - 1^2,^, -1). 

It is also clear that the thermodynamic flux j(») admits the following simple 
form : 

(2.16) j{W) = (0,p,pu)\ 

The explicit evaluation of matrix Y{v) is simple ; we have 

/I 0\ 

(2.17) Y{v) = \ -V 1 

and the orbit of the particular state W=(^p,pu,pe-\--pu^)^ Gfi is a 
parabola that is expressed by 

(2.18) G.W = \{p, p{u-v), pe ^]^p{u-vff , € IR } . 

3. Invariance 

Hypothesis 3. 

The mathematical entropy is invariant under the Gahleo group. 

(3.1) VpeG, VVFel], r]{g.W) = r]{yV) 
and this hypothesis implies in particular that 

(3.2) V?;elR, VWel], r]{y{v).W) = r]{W) 

(3.3) VVF € O, T]{R.W) = r]{W). 

Definition 5. 

Systems of conservation laws invariant for the Galileo group. 

The system of conservation laws (2.1) associated with a mathematical entropy 77 
and an entropy flux satisfying hypothesis 2 (relation (2.7)) is said to be invariant 
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for the Galileo group G of transformations if for each regular solution IR'^ D 
3 it, x) I — y W{t,x) G O of the conservation law (2.1) and each Galilean 
transformation g E G, the function V{0^ ^) defined by the relation 

(3.4) ir2 d g{E^) 3 {e, ^ v{e, ^ = {g.w)ig-\e, 0) e o 

is also a regular solution of the conservation law (2.1). It is the case in particular 
when g = is a special Galilean transform : 



(3.5) 



d 



-[Y{v).W{O,^ + v0)] + 



d 



u{Y{v) . W{0, ^ + vO) {Y{v) . W{e, ^ + vO)) 



+j{Y{v).w{e,^ + ve)) 







— r^(Y{v).W{d,^ + vd)) + 
d 



(3.6) 



u{Y{v).W{e, ^ + vO)) r]{Y{v).Wie, ^ + vO)) 







and when = g is the reflection operator 



(3.7) 



^ 'u(R.W(0, -0) R.W(0, -i) + j{R.W(e, -0) 







(3.8) ^r](^R.w{0, -o) + -^[u{R.w{e, -0) v{R^w{e, -0) 



d_ 







Proposition 3. Transformation of the velocity field. 

We have the following properties for the velocity field associated with a special 
Galilean transformation and with the space reflection when the hypotheses 1, 2 
and 3 are satisfied : 

(3.9) Vi;GlR, yWen, u{Y{v).W) = u{W) - v 

(3.10) ^Wen, u{R.W) = -u{W). 



Proof of Proposition 3. 

We first consider the elementary calculus that express the partial derivatives on 

each side of the Galilean transformation : ^ = ^+ ^^5^ = ^ 

we remark that the relation (3.2) implies 

(3.11) v{Y{v) • W{e, <e + vO)) = r](W{e, <e + vO)) . 

We develop the left hand side of relation (3.6). We get 
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dW 
~dt 



+ V 



dW 



dx 



d 



u{Y{v).w{e,i + ve)) 'q{Y{v).w{e, i + vd)) 



d 



dW 



dx 



d 



+ ^ [u{Y{v) . W{0, ^ + ve)) ri{W{0, ^ + vO))] 

-\ 

dx I 



due to (3.11) 



-u{W) + V + u{Y{v).W) r]{W{e,^ + ve)) } = 



This last expression is identically null for any regular solution W{t,x) . In 
consequence the coefficient in front of ri{W{6, ^ + v6)) is null (see e.g. Serre 
[Se82]) and this fact is exactly expressed by the relation (3.9). In order to prove 
the relation (3.10), we develop the left hand side of the relation (3.8). We obtain 



dW d 



u{R.w{e, -0) r^iw{e, -0) 
f^\u{w{e, -0) r^{w{e, -0)1 + i-,\u{R.w{e, -O) r^{w{e, -O) 



d 



d_ 



[ [u{w{e, -0) + u{R.w{e, -0)] r^{w{e, -O) } = o 



Then the bracket is null as above and the relation (3.10) is established. 



□ 



Proposition 4. Transformation of the thermodynamic flux. 

Let (2.1) a system of conservation laws satisfying the hypotheses 1 to 3 and 
invariant for the Galileo group. Then for a special Galilean transformation 
and the space reflection q we have 

(3.12) V^;€lR, VM^eO, j{Y{v).W) = Y{v).j{W) 

(3.13) VI^gO, j(R*W) +R.j{W) = 0. 

Proof of Proposition 4. 

For a special Galilean transformation, we have from the relation (3.5) : 

d 



Y{v) 



dW dW 



dt 



dx 



+ ^ { u{Y{v) . W)) {Y{v) . W) + j{Y{v) .W)] = 



-I 

dx I 



Y{v).[u{W)W + j(W")]} + — {^vY{v).W^ + 



+ ^ { u{Y{v) . W)) {Y{v) . W) + 3{y{v) . W) } 
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= ^{ [-u{W) + v + u{Y{v).W)\ {Y{v).W) + 



+ 



[-Y{v).j{W) + j{Y{v).W)\ } 



= ^{ [-Y{v).j{W) + j{Y{v).W)\ } due to (3.9) 

= according to the relation (3.5). 

Then the relation (3.12) is established ; the end of the proof is obtained by 
writing that the conservation law (3.7) is satisfied for RmW : 



u{R.W) {R.W) + JiR^W) 



R 



dW 



dt 



d_ 
dx 



-u{W) {R.W) + 3 {R.W) 



dx 

d_ 
dx 



u{W)W + j{W) 
R*j{W) + j{R.W) 



d_ 
dx 



u{W)R.W - j{R.W) 







Then the relation (3.13) is established and the proposition 4 is proven. 



□ 



4. Null- VELOCITY manifold 
Definition 6. Null-velocity manifold. 

We denote by Qq the null- velocity manifold, i.e. the set of all states W E 
whose associated velocity u{W) is equal to zero : 

(4.1) W efto if and only if u{W) = 0. 

We remark that the denonination of manifold is appropriate because the function 
W I — y u{W) G IR is regular. 

Proposition 5. Decomposition of the space of states. 
The open cone Q is decomposed under the following form 

(4.2) n = u Y{v).no 

veJR 

and the sets Y{v) • Qq have no two by two intersection : 

(4.3) yv,wen, v^w =^ {Y{v).no) n {Y{w).no) =0. 
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In other words, the cone O is a boundle space with basis Oq? fiber G«0o over 
the manifold Oq and projection 11 defined by 

(4.4) n3W^ U{W) = Y{u{W)) .W € Oo . 
Proof of proposition 5. 

The relation (4.2) is a consequence of the following remark : 

(4.5) Y{u{W)).W e Oo 

that takes into account the relation (3.9) and the definition 6 of Oq- In a similar 
way, Y{v) • Oq is also the set of all the states having a velocity exactly equal 
to —V ; then assertion (4.3) is clear. The end of the proposition consists simply 
in using the vocabulary of topologists. We refer the reader for example to the 
book of Godbillon [Go71]. □ 

Proposition 6. The null velocity manifold is of co-dimension 1. 

(4.6) yWeQo, dimT^Qo = m-1. 

Proof of Proposition 6. 

We start from the relation (3.9) : u(Y{v) • = u{W) —v and we derive this 
expression relatively to the variable v. We obtain du {Y{v) • • (^dY{v) • = 
— 1 , then we consider the particular value v = 0. It comes : 

(4.7) \/Wen, du{W).{dY{0).W) = -1. 
For W e ^ and r e IR"^'* , we have 

du{W) . { r + {du{W) . r) (dF (0) .W)} = 

= {du{W).r) + {du{W).r){-l) = 
Then the vector r e IR"^'* can be decomposed under the form 
r = p -\- [dY {0) • W) with du{W) • p = 0. If we suppose now that the state 
W belongs to the null- velocity manifold , the condition dM(VF) • p = 
implies that p G . We deduce from this point the following decomposition 

of space IR"^'* : 

(4.8) IR"^'* = T^Qo + M{dY{0).W), W e Qq . 
The decomposition (4.8) is in fact a direct sum. 

If r e (T^r^o) n (lR(dy(0) •W)) , there exists some scalar € IR such that 
r = iJ,dY{0)»W. Then, the property that r € T^fig implies du{W)»r = 0. 
Due to the relation (4.7), we deduce that p, = and the vector r is null. Then 
the property (4.6) is established. □ 
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Hypothesis 4. 

Null-velocity manifold is invariant by space reflection. 

The mill- velocity manifold is supposed to be invariant point by point by 
space reflection : 

(4.9) VVTGfio, R*W = W. 
Definition 7. Decomposition of space. 

We introduce the two eigenspaces associated with the reflection operator R : 

(4.10) = {VreH^'*, R.W = W} 

(4.11) = {W eM'^^\ R.W = -W} . 

An immediate consequence of the property (1.13) is the decomposition 

(4.12) IR"''* = A^ © A_^ . 

Proposition 7. Constraint for the thermodynamic flux. 

(4.13) yweno, j{W)eA_^. 

Proof of Proposition 7. 

It is an immediate consequence of the relation (3.13) : j{R»W) + j{W) = 
and of the hypothesis 4 : R»W = W when € Oq . □ 

Remark 1. Linear geometry. 

The hypothesis 4 can also be written as 

(4.14) Oo C A^ 

and the null- velocity manifold is flat. 

Example. Gas dynamics (ii). 

In the case of the Euler equations of gas dynamics, we have 

(4.15) Ai = {(/?,0,e)Sp€lR, e€lR} 

(4.16) A_i = { (0, (J, 0)S (J e H} . 

We remark also that all the matrices Y{v) have a common eigenvector that 
generates a linear space Fi of dimension 1 included in Ai : 

(4.17) Ti = {(0, 0, e)\ e e IR} C Ai. 

Moreover, the half manifold = { (0 , , e)^ , e > } is a part of the boundary 
of Oq which is composed by (unphysical ?) states without any matter (p = 0), 

undeflned velocity = — = full of energy (e > 0) ! 
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5. The case m = 1 

Proposition 8. 

It does not exist any hyperbolic equation (m =1) invariant for the Galileo group. 
Proof of Proposition 8. 

The operator R is linear IR — y JR and satisfies = Id. If dimA-j^ = 1, 
then u{RW)+ -\-Ru{W) = implies that u{W) = for all W e Q. Then 
T^fto = IR and this property is in contradiction with the proposition 6. If 
dimA_-|^ = 1, i.e. RW = —W, we have from (1.12) and the derivability of the 
function IR 9 i — y Yiv) : Y{v) = a'" for some a > 0. On the other hand, 
from relation (1.14) : Y{v)RY{v) = R we deduce that V^; € IR, Y{v)'^ = 1. 
Then a = and Y{v) = 1. But this property is in contradiction with the 
property (3.9) : u{Y{v) W) = u{W) — v . The proposition is established. □ 



6. Galilean invariance for systems of two conservation laws 
Theorem 1. 

When Q c ]0, +oo[ x* IR C IR^'* , a system of two conservation laws invariant 
for the Galileo group is parameterized by the scalars a > , /3 > and by a 
derivable strictly convex function ]0, +oo[ 3 ^ i — > cr(^) € IR. It takes one of 
the following forms : 

(i) Hyperbolic Galileo. 

We have for this first case 

(6.1) a'iO < 

and the space of states O is included in the following one : 

(6.2) = = (^^^ gir2'*, e>o, ICI<y|l^l}. 

The system of conservation laws takes the algebraic form 




with a velocity u{») given by the relation 

and a function n(«) named here the mechanical pressure and satisfying the 
relation 
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if we denote by cr*(») the dual function of cr(»). Moreover, the function r]{») 
defined by 

(6.6) ri{e, = ^(\/^'-^) 

is a mathematical entropy associated with the hyperbolic system (6.3). 

(ii) Elliptic Galileo. 

In this second case, we have 

(6.7) a'(O>0. 

The elliptic Galileo system of conservation laws admits the expression 




with a velocity w(»), a mechanical pressure !!(•) and a mathematical entropy 
r]{») defined by the relations 

(6.11) rj{e, c) = 




Proof of Theorem 1. 

• We have = Id in the linear space H^'*. If dimA^ = 2, then R = 
Id, = {0} and due to the relation (4.13), j{Wo) belongs to A_^ when 
Wo lies in Qq. Then j{Wo) = if Wq e Qq- We deduce from (3.12) that 
j(Y{v)*Wo) = Y{v)*j{Wo) for each Wq e Qq- Then, according to (4.2) 
and the preceding point, we have j{W) = for each W G Q. We deduce 
from the proposition 2 and the property (2.9) that du = and this fact 
contradicts the relation (4.7). Then dimA^ < 1. Moreover the unidimensional 
(due to (4.6)) flat manifold T^Qq is included in A^. Then dimA^ > 1 and 
dimA^ = dimA_^ = 1. 

• We differentiate the relation (1.14) relatively to the variable v : 
dY{v).R.Y{v) -\- Y{v).R.dY{v) = and we take v = : 
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(6.12) dY{0).R + R.dY{0) = 0. 

For r € , we have dY{0) •r + R* dY{0) • r = and 

(6.13) VrGA^, dF(0).r€A_i. 

In a similar way, for r G we have — dy(0)»r + R»dY{0)*r = and 

this imphes 

(6.14) VrGA_^, dY{0).reA.. 



• Let (r+, r_) be a basis of the hnear space IR ' composed by a non null 
vector r+ of A^ and a non null vector r_ of A_^. We define a>0 by the 
condition dY{0) -r^ = —ar- after an eventual change of the sign of r_. Due 
to (6.14), the vector dy(0)«r_ belongs to the linear space A^ and can be 
written under the form : dF (0) • r_ = /3r^ . If the scalar a is null, we can 
express the matrix Y{v) as 

^0 



Y{v) = e^d^(O) = exp v 



(6.15) -I yiy y ^^^]^ \ 

and for Wq = [Oq, 0)* G f^o , we have 
relation (4.2), this property implies that 




Y{v).Wo = {00, Oy. Due to the 
Q is included inside the subspace 
A^, that contradicts the definition 2 that claims that O is an open set of 
IR2't. Then a>0. 



If the scalar (5 is null, we can express the matrix Y{v) as 



(6.16) 



Y{v) = e^d^(O) = exp v 




with 



Then for Wq 

u{Y{v) •Wq) = —V due to the relation (3.9). We deduce the expression 
u{W) = (/6 for the velocity field. Moreover due to the invariance (3.2) of 
the mathematical entropy, we have the following calculus : 

77(^0, -^^0) = r]{W) = r]{Y{v).Wo) = r]{Wo) = 77(^0, 0) 

and the mathematical entropy is function of the unique variable 6. In conse- 
quence the function •) can not be a strictly convex function of the pair 
(^, Q. Due to the general choices done in the section 2, this case must be ex- 
cluded and the matrix of the operator dy(0) has relatively to this basis one 
among the two following expressions : 

' -p^ 

-a 



(6.17) dy(0) = 

(6.18) dF(0) = 



a> 0, /3 > 







-a 



a>0, /3>0 
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• Case (i). When the operator is defined in the basis (r+, r_) € x A_^ 
with the matrix (6.17), the end of the construction of the system of conservation 
laws can be done as follows. We first remark that 

(6.19) (dF(0))^ = a (3 Id; 

we have the expansion 

2 3 

Y{v) = e^^^(O) = Id + ^;dy(0) + ^a/3Id + ^a/3dy(0) + ••• 
and the sum of the previous series is equal to 

^^9^ V(.A - ( ch(?;v^) -vW« sh(?;v^)\ 

[p.ZK)) - l^_y^gh(i;v^) c\i(v^) )■ 

For a state Wq = (^o, 0)* G f^o , we have due to the expression (6.20), 

with u{W) = —V. We deduce the relation (6.4) with the adding condition 

ch ( u{W) y/a p ) V a 

• We focus now on the mathematical entropy. We first note (t(«) the restric- 
tion of the mathematical entropy !]{•) tothesubset Qq = f2n (]0, -|-oo[x*{0}). 
With the preceding notations, we have necessarily from the hypothesis (3.2) 
ri{W) = ri{Wo) = (t{6o) that establishes exactly the relation (6.6). A natural 
question is to verify that the mathematical entropy !]{•) is a strictly convex 
function, when cr(«) satisfies the same property. We set 



(6.22) ^ = W^2_^ 
V a 

and we have from the relation (6.6) the following calculus : 

CAC aAa ^ /-AA dr] e , drj 13 ( , 
Cd^ = Ode CdC, ;SZ = 7^' ^ = 7^ 

9^ = --^- ■ Then 

det(d^) = ^ 



56*2 \dedc 
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function r]{») is scritly convex. 

• For Wo = (6'o, 0)* G Oq , we have from (4.13) : j{Wo) G A_^ , then we can 
write it under the form : 



(6.23) VT^o = (^0,0)* gOo, j{Wo) = (0,n(^o))* e A_i , 



where the function ]0, +oo[ 3 9q i — )■ n(^o) ^ IR- remains to be determined. In 
order to find the algebraic expression of the thermodynamic flux J^*), we derive 
now the relation (3.12) relatively to the variable v. It comes : 

dj{Y{v).W).diY{v).W = dY{v).j{W) 

and taking the particular value ^; = : 

(6.24) VVTgO, dj{W).dY{id).W = dY{0).j{W). 

We apply the compatibility condition (2.9) to the vector dY{0) • Wq with Wq G 
Oo- Taking into account the relation (4.7), we get 

<f.dY{0).j{Wo) = ip.dj{Wo)*dY{0).Wo = 

= -r]*{ip)du{WQ).dY{0).Wo = r]*[if) id est 

(6.25) VVFoGl^o, ifdY{<d).j{Wo) = T^*{if) with = d77(VFo) • 
We have also from the relations (6.6) and (6.22) : 



and if W = Wq = (Oq, 0)* G Oq , we have (f = (cr', 0). We have also the 
following calculus : 



We have also 'r/*(a'(0, 0) = - 7]{^) = - ct{^) = a*(a'(0) and 



due to the relation (6.25) and the preceding development, we have necessarily 






(6.27) n(^o) 
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in coherence with the relation (6.5). Then the thermodynamic flux function 
j{W) can be easily deduced from the relation (3.12). We get, due to the condi- 
tion V — —u{W) : 

/ ch (u ^/a /3) 



j{Y{v).Wo) 



\fW[oL sh {u v/q; /3) 
\/oc/ j3 sh (m ch{u yja j3 ) 



= n(^o)ch(.v^) (VN^^C/^ 



n 



^2 



PC/a 



.n(^o)) 

due to (6.4) 
thanks to (6.21) 



and the expression (6.3) of the hyperbolic system is established in this first case. 

• We still have to verify the global coherence of what have been done, i.e. 
that the function ri{») introduced at the relation (6.6) is really a mathematical 
entropy for the system (6.3). We first have th.{u^/Q~P) = \f^JoL C,/9 then by 



derivation 



y/a j3 du = J ^ (^9 dC, — C, dO^ and due to the relation 



ch2(nA/a/3) 



(6.21) we have the following expression for the derivative of velocity : 
(6.28) adu = (^dC - Cd^) • 

Let Wit., x) = {6{t, x), ({t, x)y be a regular solution of the system (6.3). We 
have, with the notation (6.22) : 



dr] d 
dt dx 



d 



du 
dx 



da / d d 
d^ \dt dx 



du 
dx 



9 fd9 ^ d9^ 13 C (dC dC\ 

a ^ \ dt dxJ 



du 
^ dx 



a 
1 



= 9 



dx dx\^ a 



a 



^du ^ _^/n^ 

dx dx\^ 



^ V dx a ^ 



C--9^ 
dx dx 



du 



du 

^ dx 
due to (6.3) 

due to (6.22) 



(-4-'(0 + ^)|-/3.'(0 n(og 



due to (6.28) 
due to (6.5) and (6.6) 
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= 

by definition of the dual of a function in the sense of Moreau. Then any regular 
solution of the system (6.3) is also solution of the equation (2.2) of conservation 
of the entropy ; in other terms, the relation (2.9) is identically satisfied. This 
fact ends the first part of the theorem 1 where has been developed the case of a 
jacobian matrix dy(0) given by the relation (6.17). 

• Case (ii). When the matrix dy(0) is given by the relation (6.18), we 

have 

(6.29) (dy(0))^ = -a 13 Id. 
Then 

Y{v) = e^dr(O) 

2 3 4 

= Id + vdY{0) - ||-a/3Id - ^a/3dF(0) + ^{^Pf Id + • • • 



= cos(v J~a^ ) Id H sin(7; J~a^ ) dy(0) 



(6.30) Y{v) = 



cos {y \Ja^ ) yW/oi sin {y \Ja (3 ) 



— \/ Oij f3 sin {y \Ja^ ) cos {y \/al3 ) J 
For Wo = (6*0, 0)* G Oq , we have 

We deduce that necessarily the relation (6.9) is compatible with the previous 
expression and 



(6.31) 00 = . = V^' + ^ if ^ = 0* €0. 

cos [u{W) \Ja p ) V a 

• As in the hyperbolic case, the condition u{W) = implies that the 
second component C of the state W is null. Then the necessary condition 
r]{W) = (j{6o) conducts to fix the entropy on the null velocity manifold Qq as 
a strictly convex function of some variable > ; we set : 



(6.32) r]{e, = ^(0 , with i = \ d^ + ^ . 

V a 

We have again to determine the condition for the convex function (t(») to 
construct a stricly convex entropy rj from the relation (6.32). We have : 

^d^ = ^d6' + -CdC, ^ = 7^' 7^ = ~c^' 
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W^ae" +^^- • Then 

when the condition > of relation (6.7) is true. In these conditions, we 

have also > , and we have established that the function r]{») is a scritly 

convex function of the pair [6, C)- 

• We search now the thermodynamic flux j{Wo) on the form (6.23). The 
determination of the entropy variables is easy : 

^ P 
J 

and ^{Wq) = {a', 0) when Wo € r^o- We have now 
^.dF(0).i(iyo) = (<.',o) (I («) = (O, («) 

= pa'(eo)u(eo) 

and the relation (6.10) is established due to (6.25) which remains true for this 
second case. The end of the determination of the thermodynamic flux j{W) 
for an arbitrary state W is easy, taking into account the relation (3.12) and the 
expression (6.30) : 

j{W) = Y{-u{W)).j{Y{u{W)).W) = Y{-u{W)).j{Wo) 

/ cos{u ^/a f3 ) —y//3/a sm[u^/aP)\ f \ 



(6.33) WW = {e^cYen, ^ = J (^e, ^c) 



\ ^/a/~P sin(^u cos(w ^/a P ) J \Il{6o) J 

= n(flo)cos(«^) (-vW^vW^C/^^ 

= n('./fl2 _|_ fl^2/^ ] ^ I ^f'^" I thanks to (6.31) 

V ^ y^iT^z V J 

and the algebraic expression (6.8) of the hyperbolic system is established. 
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• As in the first case, we must verify that the candidate !]{•) for a mathemat- 
ical entropy is in fact a correct one, i.e. that the quantity ri{W) is advected 
with the velocity u{W) for the regular solutions of the equation (6.8), with 
the entropy flux u{W) ri{W). Taking into account the relation (6.9), wc have 

tg{u^aP) = ^J^ja C/e then 



and due to the relation (6.31) we have the following expression for the derivative 
of velocity : 

Let W{t, x) = {6{t, x), C(t, x)Y be a regular solution of system (6.8). We 
have, with the notation (6.32) : 

uT) / d \ Ou 

i + rM^ = (a* + "^"^ + "9^ 

da / d d \ du 



^ \dt dxJ a ^ \ dt dx 



+ 7] 



du 
dx 



a 
1 



= - e 



ox ox \t, a 



a 



dx dx\^ 



C V dx a ^ 



dx dx 



du 



\ du 

due to (6.8) 
due to (6.32) 



= (-^.'(0 + ^(0 + -* HO)) £ 

= 0. 

The proof of Theorem 1 is completed. 



due to (6.34) 
due to (6.10) and (6.11) 

□ 



Remark 2. About the system. 

Following a remark proposed by F. Coquel in march 2000, we can fix /3 > 
and take the limit a — > +00 for the systems found at the theorem 1. Then 
the velocities u{W) defined in (6.4) and (6.9) tend to zero and both systems of 
conservation laws (6.3) and (6.8) admit the following formal limit 





n(^) 



0; 
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wc observe here that this hmit system is not the p-system that takes the classical 
form (see e.g. [GR96]) : 

«' l(0*sU) = » 

Remark 3. Preliminaries. 

In all fairness, the elliptic Galileo version of two by two Galileo group pre- 
serving systems of conservation laws satisfies the relation (3.9) u{Y{v) •W^ = 
= u{W) — V only for sufficiently small velocities v. It is clear from the relation 
(6.9) that the velocity u{W) should be defined modulo some additive constant 
(equal to ir/^/a /3 ) because only the expression tgi^^/ a /3 u{W)) is well defined. 
Then the hypothesis 2 should be adapted to systems of conservation laws whose 
velocity belongs to some quotient group of the type 'iR/dj.TL). These devel- 
opments have not been realized at this moment [november 2000], and this fact 
explains the word "preliminary" in the title of our contribution. This remark is 
also to be done for the three by three "elliptic Galileo" system developed in the 
next section. 



7. Galilean invariance for systems of three conservation laws 
Theorem 2. 

When O C ]0, -|-oo[ x* IR^'* C IR^'* , a system of three conservation laws 
invariant for the Galileo group has one of the three following types : hyperbolic 
Galileo, elliptic Galileo or nilpotent Galileo. The hyperbolic system is parame- 
terized by the scalars a > , 5 > and by a derivable strictly convex function 
]0, -|-oo[ X IR 3 {a, j3) I — )■ cr(a, /3) € IR; we denote by cr*(«) the dual function 
of a{.). 

(i) Hyperbolic Galileo. 

We have for this first case 

(7.1) Va>0, V/3GIR, — (a, /3)<0 

da 

and the space of states O is included in the cone defined by 

(7.2) 0+ = = (^, C,V^)*€IR^'\ ^>0, Kl<yfl^l}- 
The system of conservation laws takes the algebraic form 
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with a velocity u{») given by the relation 



(7.4) u{W) = 



'06 



argth I 



a e 



and a so-called mechanical pressure function !!(•) satisfying 

1 a* {A, B) 

~b A ' 
Moreover, the function ri{») defined by 



(7.5) n(^, /3) 



with {A, B) = da{^, /3) 



(7.6) 7/(^,C,^) = ^{\0'- 



a 



is a mathematical entropy associated with the hyperbolic system (7.3). 

(ii) Elliptic Galileo. 

We have in this second case 

(7.7) Va>0, V/3GIR, ^{a,P)>0. 

The elliptic Galileo system of conservation laws admits the expression 



d 







d 



(7-8) ^ C \ +^{u{W) C + 



dt 



dx 



with a velocity w(»), a function n(») and a mathematical entropy ?7(«) 
fined by the relations 



(7.9) 
(7.10) 



u{W) = 



arctg I 



ah 

1 a* {A, B) 
b A 



a e 



(7.11) 77(^, C, V^) = ^(^V^' + 

(iii) Nilpotent Galileo. 

For this third case, we suppose 



bC 



a 



with {A, B) = da(^, 

A- 



|(«,/^)<o 



(7.12) Va>0, V/5G1R 
The nilpotent Galileo system of conservation laws takes the algebraic form 

/e\ . , //)\ _ „ / \ 



(c)^£H(c).n(..-A5) 



1 

\a6 J 
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with a velocity field u{») given by the relation 

(7.14) u{W) = j- 

and a mechanical pressure !!(•) satisfying 

(7.15) n(a, (5) = -i , with (A, B) = da{a, (5) 

B 

or in an equivalent way 

1 ^*(^("' ^)) 

(7.16) n(a. fS) = g— — . 

Moreover, the function r]{») defined by 

(7.17) v{e,c^) = ^(^'V^-^y) 

is a mathematical entropy associated with the hyperbolic system (7.13). 
Remark 4. Pressure and duality. 

The case of gas dynamics corresponds to 6 = p, if) = pE^ a = h = 1 ior 
the nilpotent Galileo system of conservation laws. The entropy (t(«) at null 
velocity is defined from the classical thermostatic entropy function E(«) which 
is concave and homogeneous of degree one relatively to the extensive physical 
variables of mass M , volume V and internal energy 8 [Du90] . Introducing the 
intensive thermostatic variables of temperature T, thermodynamic pressure p 
and massic chemical potential we have the classical fundamental relation of 
thermostatics (see e.g. Callen [Ca85]), 

(7.18) = T dE(M, V, 8) - pdV + fidM 

and taking into account the Euler relation for homogeneity of degree one 

(7.19) 8 = TT. - pV + pM , 
we have necessarily : 

(7.20) 77(p,0,^) = -E(p,l,^) = a(p, ^). 

We take V = 1, M = p and 8 = ip\ we deduce from the relations (7.19) and 

(7.20) the identity : 

(7.21) a = iv^- I 

and by application of the relation (7.18) in the particular case V = 1 , we 
deduce : 
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(7.23) <-*(f,-^) 



(7.22) da = ^dp- ^d'^. 

Both identities (7.21) and (7.22) establish that we have 

f 1^ 1 \ P_ 
T ' 

In this context, the relation (7.12) can be written as ^ > and the following 
calculus, issued from the relation (7.16) : 

m /Ml 
np,i') = -JZYJT) = (IM = At'-t 

gives an intrinsic definition of the mechanical pressure n(») as identical to the 
thermodynamic pressure p{») with the help of the relation (7.23) in terms of the 
dual of the thermostatic specific entropy. We remark also that, in some sense, 
the theorem 2 establishes theoretically the Euler equations of gas dynamics. 

Proof of Theorem 2. 

• We derive the relation (1.14) relatively to the velocity v : 
dY{v).R.Y{v) +Y{v).R.dY{v) = 

and we consider the particular case = 0. It comes : 

(7.24) dy(0).i? + i?.dr(0) = 0. 

Then the range of the linear space by the operator dy(0) is included inside 
the eigenspace A_^ and the range of the linear space A_j by the same mapping 
is included inside the eigenspace A_^. Moreover, the kernel of dy(0) is stable 
under the action of the symmetry operator R : 

(7.25) dy(0)(AJ c A_^ 

(7.26) dy(0)(A_J c A^ 

(7.27) R(kerdY{0)) C kerdy(O). 

We deduce from the relation (4.6) (dim T^Q,q = 2 in our case) and from the 
relation (4.14) {Qq C A^) that we have necessarily : 

(7.28) dimA^ > 2, T^Oq C A^ . 

• The case where dimA^ = 3 is not possible. Indeed we would have 
A_^ = {0} and the relation (4.13) establishes that in this case the thermo- 
dynamic flux j(») is identically null on the manifold flo then on the entire 
space ft. Moreover, the relation (3.10) joined with the relation R = Id shows 
that u{W) = for each W G Q. This fact contradicts the relation (4.6) that 
claims that dim T^Qq = 2. Taking into account the relation (7.28), we have 
established that we have necessarily 
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(7.29) dimA^ = 2, dimA_^ = 1. 

If dimkerdy(O) = 3, then Y{v) = Id for each real parameter v and the 
Gahleo group does not operate anymore on the cone O. In particular, we have 
a contradiction with the property (3.9) : u{Y{v) •) W) = u{W) — v. 

• If dim kerdy(O) = 2, due to the stability property (7.27), we are neces- 
sarily in one of the two following cases : 

(7.30) kerdF(O) = 

(7.31) 3 ro e A^ , ro 7^ 0, 3 r_ e A^ , r_ 7^ 0, ker dY{0) = span < ro , r_ > . 

If the relation (7.30) holds, there exists a basis (r+ , ro) of the linear subspace 
A-j^ with dy(0)«r+ = dy(0)«ro = 0, a non null vector r_ e A_-|^ and a 
scalar a > such that inside the basis (r+ , r_ , ro) , the operator dy(0) 
admits the following expression : 

/O -a 0\ 

(7.32) dy(0) = in the basis (r+ , r_ , ro) . 

\0 0/ 



that belongs always in the null manifold Oq and we have a contradiction exactly 
as in the preceding point. If the relation (7.31) is active, there exists a scalar 
a > such that inside a basis (r+ , r_ , ro) consructed as previously, we have, 
taking into account the relation (7.25) : 




and for Wq = 




we have Y{v) • Wq 





/ 
-a 
\ 




in the basis (r+ , r_ , ro) . 



Then 




and for 
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We deduce from the relation (3.9) that that u{6,(,iIj) = (/{a9) and for W = 
(^, i/jY G r^, we have necessarily ri{6, ''P) = cr(^5 i^) where a{») is the 
restriction of the entropy !]{•) to the null manifold D.q. Then the mathematical 
entropy ri{») cannot be a strictly convex function and this case has to be 
excluded. In consequence we have necessarily dim (ker dy(0)) < 1 . 

• We observe now that we have also necessarily dimkerdy(O) > 1 because 
taking into account the relation (7.25), the operator dY{0) can be considered 
as a linear mapping between the linear space of dimension 2 and the linear 
space A_^ of dimension 1, therefore with a necessarily non degerated kernel. 
According to the previous point, we are inconsequence in the particular case : 

(7.34) dim (ker dF(0)) = 1 , ker dF(0) C A^ . 

With analogous notations as in the previous subsection, we introduce a non null 
vector ro {vq G A^ ) that generates the linear space ker dF (0) . We complete 
the basis of A^ by some non null vector r_|_ and due to the inclusion (7.25), its 
image r_ by dy(0) is necessarily a non null vector of the eigenspace A_^ : 

(7.35) 3a >0, dF(0).r+ = -ar_ . 

We write now the vector dY{0) •r- € A^ inside the basis (r+ , ro) : 

(7.36) dF(0).r_ = cr+ - dro- 

If c ^ 0, we set r+ = — dro/c and we have dy(0)«?+ = —av- ; 
dy(0) •T- = cr+ . Inside the basis (r+ , r_ , ro) the matrix of the operator 
dy(0) takes the form : 

/ c 0\ 

(7.37) dF(0) = -a , a>0, c^^O. 

\ 0/ 

We change our notations and replace the vector r+ initially introduced by the 
new vector r+ . When the scalar c introduced at the relation (7.36) is null, we 
have necessarily d ^ and after an eventual change of the sign of the vector 
ro , the matrix of the operator dy(0) inside the basis (r+ , r_ , ro) can be 
expressed as 

/ 0\ 

(7.38) 3a > 0, 6 > 0, dy(0) = -a in the basis (r+ , r_ , ro) . 

\ -6 0/ 

• Case (i). We develop now the two particular cases (7.37) and (7.38). If 
the relation (7.37) holds, we have in a first opportunity : 

/ -6 



(7.39) 3a > 0, 6 > 0, dy(0) 



—a I in the basis (r+ , r_ , ro) 

V 
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and we are exactly in the case (i) of an "hyperbolic Galileo" system of con- 
servation laws. The proof follows what have been done at the theorem 1. By 
exponentiation of the relation (7.39), we have : 

/ ch(^v Vab) —■sjhja sh(i;-\/a6) 

(7.40) Y{y) = -Va7^sh(7; V^) ch(?;A6") 

V 

For an arbitrary state W = {0, ipY E Q we have : 



Y{v) 




ch(yVab ) — C y^b/a sh(vVab^ \ 
— y^ajb 6 sh(v y/ab^ + C ch.(v y/ab ) 



Because u{Y{u{W))*W) = and C = IR {0} x*IR, we deduce 
from the previous equality the expression (7.4) of the velocity field. We then 
have 9 ch (v \fah ) — C \fbj~a sh {y Va/b ) = ^Q'^-hCja and the necessary 
condition T\iyV^ = 'ri(Y{u{W)) ^W) shows that the relation (7.6) holds. 

• We must look now to the precise conditions that makes the function 7]{») 
defined in (7.6) a strictly convex function when the property is satisfied for the 
two- variables function a (a, (3). We set as in the proof of the theorem 1 : 



(7.41) ^ = 

and we have 

= ede 




bC 



a 



-CdC, 
a 



d( dip 
Then 



-P ,3 



da 



+ 



drj 
89 



b 

P = - 
a 



9 da 

1 da ' 



drj 



bCda 
a ^ da ^ 



drj 

di/j 



da 

d/3 



d9 dC 
9^ da 



9_Cda__ 6>C d'^a 
^ da ^ da^ 



e d^a 



a 



da 
d'^c 



e da^ 



-Pt 



(7.42) d^T/ = 



i da d/3 
( e da 

9C da 
^ da 



d'^T] 



^ e da^ 



^2 d'^a 



e da^ 
9C d^a 



9_Cdo^_ _ 
^ da ^ da^ 



da^ 



d^a 



?1 ^ 
^3 d^ 



9C d'^a 



^2 da^ 



9 d^a \ 

C da dp 

C d'^a 
-Pt 



\ 



9 d^a 
J dad 15 



-Pt 



C d'a 
i dadp 



^ da dp 
d^a 
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The two by two sub-matrix at the top and the left of the matrix of the expression 
(7.42) admits the following determinant : 



det 



l<i<2, l<j<2 



da do? 
b 1 da d'^c 



a ^ da do? 

which is necessarily strictly positive in order to get the strict convexity of the 



function r\. We know that the second derivative 



is strictly positive and 



that a b and ^ have the same property. Then we have — > and the 



relation (7.1) is established. We develop now the determinant of d'^r] by using 
the third column. It comes : 



det (^d'^rij = - 



+ - 



e d'^ 



^ dad/3 



+ 7^7 



p da d^a d'^a 
^ da do? dp"^ 

e^da 

e d'^a 



+ 



eC d'^a 



dcP- 



+ 



i dad/3 

p da d'^a d'^a 
^ da do? dl3'^ 



^ da dp 
Cda ^ 

a 



d^a 
dcP- 



e^da_ 2 
^ da ^ ^2 da^ 
_ C d'a 
^ C dad/3 
eCda eCd'^a 



e d^- 



d^a 



^ dad/3 



C 

^ dad/3 



^4 \dad/3 



^4 \dad/3 

p da d'^a d'^a 
^ da do? d/3'^ 



^ da da^ 

e 

e da ^ ^2 d'^a 

^ ^ da da^ 
d 



_p ^ ^ + ^2 ^2 ^ 
^ da da^ 

-pC 

^ da do? 

-pC, 



1 



d'^ 



da d/3 



2 da_ 

da 







pC 



^-(-p^C'+pO')+ ^ 



(p'e-p0') 



p da d'^a d'^a 
^ da da^ d/3'^ 



+ ^ 



2 da 



da d/3 
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p da 
C da 



and 

(7.43) det (d^Ty) = 



da d/3 
b 



due to (7.41) 



da 



a 



d^a 



d^ 



a 



2 1 



(7.44) if 



a y/e'^-bCya da L da^ d^^ \dad/3, 

The matrix d^r] is definite positive if it is the case for the matrix d^a and if 
1^ < 0, i.e. when the condition (7.1) is satisfied. 

• In order to obtain the algebraic expression of the thermodynamic flux j{m) , 
we first evaluate the entropy variables 

/' dr] dr] drj \ _ / 9 da bC da da \ 
Kde'dC'di^) " a^' a^' 
on the null- velocity manifold Oq that corresponds to {9 , C, , ip) = (^o , , '0o) • 

Then (pQ = { -;t—(0o, ipo) , , Tr^i^^Oi V^o) ) • Taking into account the relation 
\ da dp J 

(7.39), it comes : 

(7.45) (^o.dy(O) = (o, -6^ , o). 

We introduce the mechanical pressure n(«) as a notation : 

(7.46) V T^o = (^0, 0, ^o)* € Oo , = (0, n(^o, ^o) , O)* 

and the relation • dj + 77* dix = applied against the vector dF (0) • VFq 
gives 

6^(^o, V^o) n(^o, V^o) + o\A, B) = 

with A = ^{00, '0o) 5 B = ^{Oq, ipo), and this property is exactly the rela- 
tion (7.5). The general expression of the thermodynamic flux j{W) is obtained 
thanks to the relations (3.12) and (7.40) : 

j{W) = Y{-u{W)).j{Y{u{W)).W) 

ch[uVab ) ^/b/a sh(uVab) 

^/a/b sh(uVab^ ch (u Va 6 ) 
1 

e' 





u(^e^-be/a,^ij) 





(7.47) yw = 



C \ en, j{w) = 



^/d^-be/a 



Joined with the expression (7.4) of the velocity field, the relation (7.47) estab- 
lishes the expression (7.3) of the hyperbolic Galileo system of conservation laws. 
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• We verify now that the function !]{•) defined at the relation (7.6) is effec- 
tively a mathematical entropy associated with the flux u{W) r]{W) ; in other 
words we have the additional conservation law 
driiW) driiW) 



(7.48) 



+ 



{u{W)r]{W)) = 



dt dx 

if W{x^ t) is a regular solution of the conservation law (7.3). By differentiation 
of the relation (7.4), we have : 

(7.49) a du = {OdC-C de) . 

We have now do develop some algebraic calculus : 



dri ci ^ . /a a \ (Ju 



d 



d 



du 



da 



/ d a \^ ot, / a a \ ^ 



da L 89 
da 
da 



d 



d 



dt dx) d( \dt dx) 
9 /■ d ^ \0 ^ C f d d \ 
^ \dt dx) a ^ \dt dx) 



c 



da f d d \ , du 

d(3 \dt ^ dx) ^ dx 
due to (7.44) 



9_ da_ 
1 da 



h C da 
a ^ da 



^du d /-Q/^ ^\ ^ 
dx dx \ ' ^ 



d(5 dx ^ dx 



according to (7.3) 



da du 
da dx 



du 
dx 



da 



da du ^ ^ du 
d/3 dx ^ dx 
_ObdC_ bCd^- 
^ a dx a ^ dx. 
da 



b da , 



. dx dx. 



= 



and the property (7.48) is established. The study of the first case is over. 
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'a sm[v 
cos (v ^/ab ) 



• Case (ii). We consider again the expression (7.37) of the jacobian matrix 
dy(0) now with c > and the relation (7.37) can be rewritten as : 

/ 6 0\ 

(7.50) 3a >0, 6>0, dF(0) = -a in the basis (r+ , r_ , Tq) . 

\ 0/ 

The relation (6.30) established during the proof of the theorem 1 can be repro- 
duced without any modification and we have in consequence : 

; (v \fah ) \/b/a sin {y \fah ) 

(7.51) 



Y{v) = I -sjajh ^m.{y-Jah) 


then for (^, -0)* G O , we have 

/ ^ \ I Q cos [y \l ah ) + \Jhja ( sin (y \fab ) 
^(^)» C = -^0/6 6* sin(i; + C cos(i; 

We know that the reference basis (r_|_ , r_ , ro) belongs to the product of spaces 
X A_^ X and Oq C A^ . We deduce that the second component of the 
state Y{v)»W that belongs to the null velocity manifold is necessarily null and 
we have 

(7.52) tg(nVa6) = a/ - ^ 

and the relation (7.9) is established. We deduce naturally 



Q 

(7.53) cos{uVah) = - 



sm{u'\/ab) = 



a ^ ' 



a 



(7.54) Y{u{W)).W = U,0,ijY, e = V^' + bCVa 

and the relation (7.11) is a consequence of the invariance (3.2) of the mathemat- 
ical entropy for the transformation Y{v). 

• We consider now as given the strictly convex function (t(«) which is the 
restriction of the mathematical entropy ?7(«) to the null manifold fto. We must 
verify that the mathematical entropy is also a strictly convex function of the 
triplet {9, ip) . We set as above p = h/a and we have 

^ dr] dr] dr] ^ f da (da da \ 

Then 



dO'dC'dipJ y^da'^^da'd/SJ 



d^ 



a 



_ (^da_ _ 
^ da ^2 da^ 



d'^T] 
dOdC 



OC da OC d'^a 
^ da ^ da^ ' 
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and 



= p 



^ dadl5 



P 



^ dad/3 



(7.56) d^T] = 



da 5V 



/ c!_,_ 

^ e da e 



do?' 

eCda eCd'^a 



e d^a 



do? 



\ 



i dad/3 



^ da ^ ^2 da^ 
da 2 C 5V 

C d^a 
Pt 



e 



d'' 



a 



\ 



do? 



P 



i da dp 
C d^a 



i da dp 
d^a 



^ dad/3 



d/3'^ 



I 



This matrix is identical to the one presented at the relation (7.42), except that 
the variable p must be changed into —p. But this change of sign is exactly what 
is necessary to modify the definition of the variable ^ from (7.41) to (7.53). We 
observe that the two by two minor determinant composed by the two left lines 
and the two first columns of the right hand side of the relation (7.56) is equal to 

I 1^ 1^ and is strictly positive when |^ is strictly positive. In consequence 

the hypothesis (7.7) is clearly established. The computation of det (d^?7) done 
during the study of the hyperbolic case conducts to the relation (7.43) and the 
change of the variable "p into —j) establishes that 

(7.57) det(dr;) = --[ — —- J 

which is strictly positive when cr(«) is strictly convex and when the relation 
(7.7) is satisfied. 

• The constitution of the thermodynamic flux is obtained exactly as the 

hyperbolic case. Taking into account the relations (7.50) and (7.55), we have, 
for Wo = {00, 0, i/jy G no : 

n o / 5 

/oa 



^(^0,0, v^).dy(o) 



da\ 

\da' ' dp) 




r. , da 

'•'TP 







and if we introduce the mechanical pressure !!(•) with the relation (7.46), 
we have from the relation (2.9) : (po •dY {O) • j{Wo) + a* {—1) = 0, and 
this calculus proves exactly the relation (7.10). The general expression of the 
thermodynamic flux is given by somes lines of algebra that are consequence of 
the relations (3.12) and (7.51) : 

j{W) = Y{-u{W)).j(Y{u{W)).W) 
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\falh 



sin(u, 




cos {u^Jah ) —^Jhja sm(w\/a6) 

''ah ) cos (n \fah ) 
1 

/_ Pb rbc\ 

\ a \ a ^ 
\ / 

/ H\ 

a 



( 

V 



taking into account the relation (7.53) 



^ n(^, ^) 



V 



e 





Then the algebraic expression (7.8) of the elliptic Galileo system of conservation 
laws is established. 

• We still have to show that any regular solution of the system (7.8) satisfy 
the conservation (7.48) of the mathematical entropy, with a velocity field given 
by the relation (7.9) that satisfies in consequence : 

(7.58) adn = ^2 /^^/^ (^dC-Cd^)- 

We have very simply 
dr] , d r \ f d d 



or] C / X / C O \ du 



da 
da 
da 
da 

e da 
^ da 



de 

e / d 



d 



d 



dt "'"^ dx J 
d 



+ 



d^ 

b C / d 



dt dx) 



c 



+ 



da 
dp 



da / d 



d_ d_ 
dt dx 



d 



du 
dx 
du 



d(3 



due to (7.55) 



du 
dx 



b C da 
a ^ da 
da . , du. du 

da _^b C da ^^da 
^ da a ^ da dp 



e 



applying the equation (7.8) 



a 



+ 



d /n 



e da bC _b C da ^ 



dx \ ^ / da a a ^ da . 



b da U 



~^ a da " ^ " 



dC 

dx 



de_ 

dx 
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= taking into account the relation (7.15). 

This property estabhshes the structure of the eUiptic Gahleo system. 

• Case (iii). We detail now the third case, when the matrix dy(0) is 
given by the relation (7.38) with a > and 6 > . We first remark that the 
matrix dy(0) is nilpotent because 

/ 0\ 

(7.59) dY{Of = 0, dy(0)^ = 0. 

\ab 0/ 

This remark justifies the name "nilpotent" given for this third Galileo group 
preserving system of conservation laws. The exponentiation of the matrix Y{v) 
is easy : 

/I 0\ 

(7.60) Y{v) = -av 1 , veM, a>0, b>0. 

\ab^ -bv 1 / 

The determiation of the velocity field u{») is a direct consequence of the evalu- 
ation of the product Y{v) •]¥ and of the remark that u(Y{u{W)) • = . 
We have 

\fw = \ c \ en, Y{v) C = -ave + c 

V^/ \ab^e - bvC + iIj 

and we deduce that we have necessarily u{W) = This expression is exactly 
the relation (7.14). We deduce : 

(7.61) \/w = {e,ci^y en, y{u{w).w) = (^'O'^-^y)'- 

• The computation of the entropy variables is a simple consequence of the 
relations (7.17) and (7.61). It comes 

(7 621 <n = ( — —\ 

^' ^ \da 2ae^di3' a 9 ' 8/3 ) ' 

With the notation p = b/a, the exact expression of the Hessian of the mathe- 
matical entropy only needs some care : 

_ _ Cda 1 2C^^ d^V C da p'^ d^a 

d^ri p (2 Q'^cr (f'rj p da 2 d'^a 
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-p 



and with the complementary notation q = C^jd we get 
\ = 



(7.63) < 



da^ ^ 9 
pq da 

1 



da 1 2 4 
1 2 '^d'^o- 



pq da 

Top 

p da 

~e dp 



1 2 r,d'^a 



pq 



I 2 2 



\ 



2^ 



,d^\ 
d(3^ 

d'^a 



-pq 
^2 



5^2 
(J 



a^2 / 



a^2 ^ a^2 

The two by two minor determinant at the bottom and the right of d2?7 is equal 

to — ^f^f^gl, which establishes the necessary condition (7.12) when we have 

supposed that the domain Q is composed by states W — {9, ipY that satisfy 
the condition 

(7.64) ^ > . 

The evaluation of the determinant of the relation (7.63) is easy. We multiply 
the last line hy pq and we add it to the second line. We obtain : 

det(d277) = 

pq da 1 



d'^a 
do?- 



d'^a 



g2 1 2 4 

+ 4^ ^ W 



pq da 
~d d^ 



1 . 



d'a 

a/32 



p da 

~edp 

d'^a 



d^ 

dl3'^ 



d^ 



d(3^ 

d^a 



and we 



then we multiply the third column of the previous expression by ^pq^^ 
substract it from the first column ; in an analogous way, we multiply the last 
column hy pq and we add it to the second column of the determinant. We find 

"2 c^cr pq da 1 



det {d'^T]) 



d^a 



d^a 



q_ ^ ^ 2 

5^2 P e di3 e dp 2^^ dp'^ 



pq da 

T dp 





p da 
9 'dp 







d'^a 

dp^ 



d^a 

dp^ 



da d^a 



dp da^ 



The sign of det (d277) is the one of — when the condition (7.64) is satisfied 
then is coherent with the inequality (7.12). Then the function irj{») defined 
with the relation (7.17) is strictly convex. 



35 



FRANgois Dubois 



• The thermodynamic flux is constructed exactly as in the two previous 

cases. We use the definition (7.46) of the function !!(•) ; then the expression 
(7.62) of the entropy variables shows : 

/O 



^.dY{0).jiWo) 



/da da \ 

va^' ' op) 







-a 










0-6 

The relation (7.16) is then a consequence of the remark that dj{W) • dy(0) • W 
= dy(0) •jiW) joined with the relation (4.7). The general expression for the 
function j{W) is a consequence of the relation (3.12) and the expresison (7.61) 
of a state inside the null velocity manifold : 

2a e 
( 




\ 







2a e 



J 



^ 2a e 



due to the relation (7.14). 



The algebraic expression (7.13) of the nilpotent Galileo system of conservation 
laws is an immediate consequence of what have been done at the previous line. 

• As in the two preceding cases, we verify that the candidate (7.17) for beeing 
a mathematical entropy satisfies the relation (7.48) for regular solutions of the 
conservation law (7.13). We have 

dr] d , X / d d \ du 



or] C / X / c o \ ou 



/da 



+ 



ah o da 



u 



+ 



da 

dp 



dp) 
d 



( 



d_ d_ 

dt dxJ 



]e - b 



u 



da 
dp 



d d \ 
dt dxJ 



c 



du 
dx 

da \ f ^du\ 



da 



due to the expression (7.62) 



du dH 



^ aa / ^ du *^ ( n 6) ^ + 
dp \ dx dx ^ 'J dx 

taking into account the equation (7.13) 



d 



dp V dx dx 
du 
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du 
dx 
du 
dx 



e 



da 
da 



+ h 



2a e 



da 
d^ 



a 



bU 



da du 
d(3 dx 



/da da 



\da^ d(3 J 



= 



and the theorem 2 is proven. 



due to (7.14) 



Ail) 

2a e J 



□ 



(8.1) a{0, V^) = -e log 



0>O, 2p>0 



8. The Cemracs System 

This section describes the funny hyperbohc system of conservation laws that 
we have derived during the Cemracs 99. It corresponds to the elliptic Galileo 
system of conservation laws with a = b = 1 and the thermodynamics of the 
polytropic perfect gas. 

• Let 7>1 be a fixed real number, e.g. 7 = 7/5 for the air (di-atomic gas) 
at usual conditions of tempeature and pressure, see [Ca85]. We set 

and as considered in the remark 4, we can introduce the density p and the 
internal energy e according to 

(8.2) e = il; = pe. 

We have da = Q{d, ip) dO + x(6', ip) dip with 

(8.3) Q(^.V>) = 7+^ 

(8.4) x((l,f) = -| >0. 

The determination of the dual function of the entropy is easy. We have to solve, 
for a given pair (A, B) , the system 



7 + 



d 



(8.5) 

Then, following the relation (7.10), we evaluate the "mechanical pressure" n(«) 

n(^, ^) = ^114^ = ]L (eA + i^B- a{e, ^)) 
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= e + 
(8.6) n(^, v^) 



i-e -a) = 



(7-1)^ 



_ (7-1)^ 

(7-l)V^ 



7 - log 

Taking into account the relation (7.9), we set 

c 



(8.7) u{W) = arctg | , W = {6, il^f e Q 
and the Cemracs system is a simple re-writing of the relation (7.8) : 

/-c 



(8.8) <^ 




d 



e 



e 
\ 



w 




= 



Proposition 9. The Cemracs system is hyperbolic. 

Under the hypotheses 

(8.9) e>o, v^>o, u{^/e^Tc^,fP)>o, 

the system (8.6) (8.7) (8.8) of conservation laws is hyperbolic. 
Proof of proposition 9. 

• It is some kind of exercice for the studients. We do a direct proof of this 
proposition to improve the preceding assertions. The simplest way to do it is to 
consider the following new set of variables : 

(8.10) V ^ [e,^^^,^ 

We have 

(8.11) 
and 



n(^ Vi + e, ^) 



(8.12) 



Then 



{i+e)du = d^ = d(^) = -^d^ + idc 

(7-1)^ 



n(^ v^TTF, ^) 



7 



de 



(7-1)^ 



log 

d^ 



(7-1)^ 



^7(1+^2)7/2 

7 
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(8.13) = A.d^+ 



e J e 7-1 e 1 + ^2 ^ (7-i)6'v^' 

• We write again the two first equations of the system (8.8) : 

(8 15) + _|_ ^ _l_ Q 

5t 5a: dx dx 
We multiply the equation (8.14) by — C/^^ = —^/^ and the equation (8.15) 
by 1/6 . Thanks to the relation (8.11), we have the following evolution in time 
of the variable ^ : 

The equation for the variable (p is composed from the relations (8.14), (8.16), 
the third equation of the system (8.8) that is : 

dtb dxh du 

and from the relation (8.13). We set 
(8.18) 2/ = ^ ; 

we multiply the equation (8.14) by the coefficient (1 ^ — j — , the equation 

7 (^2 ^ (p2 

(8.16) by — —r- , the equation (8.17) by , ^ , and we add 

^ ^ ^ 7-1 6* (1 + ^2)' ^ V ^ (7_i)6>^ 

these three equations. The — term is absent, the term associated to — is 

ox ox 

equal to 



e 



and the coefficient of the term relatively to is : 

ox 

-^yU — -^y) — -^iy^ = -iy- 

V 7—1/ 7— 1 

We deduce from this calculus that with the variables V introduced in (8.10), 
the system (8.8) takes the form 
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(8.19) 
with 



+ + B{V) 



dt 



(8.20) B{V) = 



dx 





dx 







e 



e 



y 



(2 + ^2/' + 



7 



7 



y^ 



e 

-iy 



\ 



/ 



• We see in an clear way that the real number zero is an eigenvalue of the 
matrix B(y) of the relation (8.20). The two other eigenvalues A satisfy the 
equation : 



(8.21) + 



-y + 2r 



7 



7 



2/" 







The equation (8.21) has two opposite real eigenvalues when y > Q as proposed 
in our study [Du99]. We have effectively in this case : 



-y + 22/^ 



' y^ 



-y{^-yf 



7-1 

and the proposition 9 is established. 



7 



y^ < 



□ 



Remark 5. Curious Physics. 

The condition y > i.e. ip/9 > or in an equivalent way 11 > because 
^ > by convention. It corresponds to 

(8.22) n = - > 

because as claimed in the relation (7.22) : da = ^ d9 — ^ di/j. Then 
n = T — Jl mechanical pressure n(«) is the quotient of the 

physical pressure p{») divided by the chemical potential /j, which is not nat- 
ural. It is easier to interpret the Cemracs system (8.8) as an elliptic Galieo 
system of conservation laws whose associated thermostatics is obtained by ex- 
changing the roles of the variables 6 and tp = pe inside the relation (8.2). 
This corresponds to : 

(8.23) = pe, ^ = p. 

Then the relation da = — — 9 + ^ d^ and the condition (7.7) can be written 

(8.24) %--k>^ 
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corresponding to a thermostatics system with a negative temperature and asso- 
ciated with a thermostatic pressure p equal to the mechanical pressure 11 and 
given by the formula 

(7- l)pe 



(8.25) n(pe, p) = 



7 - log 



(7-l)p 
■ (pep 
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